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Escherichia coli (E. coli) bacteria govern their trajectories by switching between 
running and tumbling modes as a function of the nutrient concentration they experi- 
enced in the past. At short time one observes a drift of the bacterial population, while 
at long time one observes accumulation in high-nutrient regions. Recent work has 
viewed chemotaxis as a compromise between drift toward favorable regions and accu- 
mulation in favorable regions. A number of earlier studies assume that a bacterium 
resets its memory at tumbles - a fact not borne out by experiment - and make use 
of approximate coarse-grained descriptions. Here, we revisit the problem of chemo- 
taxis without resorting to any memory resets. We find that when bacteria respond 
to the environment in a non- adaptive manner, chemotaxis is generally dominated by 
diffusion, whereas when bacteria respond in an adaptive manner, chemotaxis is dom- 
inated by a bias in the motion. In the adaptive case, favorable drift occurs together 
with favorable accumulation. We derive our results from detailed simulations and 
a variety of analytical arguments. In particular, we introduce a new coarse-grained 
description of chemotaxis as biased diffusion, and we discuss the way it departs from 
older coarse-grained descriptions. 



Author Summary 

The chemotaxis of Escherichia coli is a prototypical model of navigational strategy. The 
bacterium maneuvers by switching between near-straight motion, termed runs, and tum- 
bles which reorient its direction. To reach regions of high nutrient concentration, the run- 
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durations are modulated according to the nutrient concentration experienced in recent past. 
This navigational strategy is quite general, in that the mathematical description of these 
modulations also accounts for the active motility of C. elegans and for thermotaxis in Es- 
cherichia coli. Recent studies have pointed to a possible incompatibility between reaching 
regions of high nutrient concentration quickly and staying there at long times. We use 
numerical investigations and analytical arguments to reexamine navigational strategy in 
bacteria. We show that, by accounting properly for the full memory of the bacterium, this 
paradox is resolved. Our work clarifies the mechanism that underlies chemotaxis and indi- 
cates that chemotactic navigation in wild-type bacteria is controlled by drift while in some 
mutant bacteria it is controlled by a modulation of the diffusion. We also propose a new 
set of effective, large-scale equations which describe bacterial chemotactic navigation. Our 
description is significantly different from previous ones, as it results from a conceptually 
different coarse-graining procedure. 



Introduction 



The bacterium E. coli moves by switching between two types of motions, termed 'run' 
and 'tumble' [1]. Each results from a distinct movement of the flagella. During a run, flagella 
motors rotate counter-clockwise (when looking at the bacteria from the back), inducing an 
almost constant forward velocity of about 20 fim/s, along a near-straight line. In an environ- 
ment with uniform nutrient concentration, run durations are distributed exponentially with 
a mean value of about tr = Is [2]. When motors turn clockwise, the bacterium undergoes 
a tumble, during which, to a good approximation, it does not translate but instead changes 
its direction randomly. In a uniform nutrient-concentration profile, the tumble duration is 
also distributed exponentially but with a much shorter mean value of about tt = 0.1s jsj]. 

When the nutrient (or, more generally, chemoattractant) concentration varies in space, 
bacteria tend to accumulate in regions of high concentration (4, 5|. This is achieved through 
a modulation of the run durations. The biochemical pathway that controls flagella dynamics 
is well understood [gUs] and the stochastic 'algorithm' which governs the behavior of a 
single motor is experimentally measured. The latter is routinely used as a model for the 
motion of a bacteria with many motors [l| 9^12 ] . This algorithm represents the motion 
of the bacterium as a non-Markovian random walker whose stochastic run durations are 
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modulated via a memory kernel, shown in Fig. [TJ Loosely speaking, the kernel compares 
the nutrient concentration experienced in the recent past with that experienced in the more 
distant past. If the difference is positive, the run duration is extended; if it is negative, the 
run duration is shortened. 
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FIG. 1: Bilobe response function of wild-type E. coli used in the numerics in Fig. [3l For the sake 
of computational simplicity, we have used a discrete sampling of the experimental data presented 
in Ref. [19] instead of working with the complete data set. This did not affect our conclusions. 




In a complex medium bacterial navigation involves further complications; for example, 
interactions among the bacteria, and degradations or other dynamical variations in the 
chemical environment. These often give rise to interesting collective behavior such as pattern 



formation 
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14j . However, in an attempt to understand collective behavior, it is imperative 



to first have at hand a clear picture of the behavior of a single bacterium in an inhomogeneous 
chemical environment. We are concerned with this narrower question in the present work. 

Recent theoretical studies of single-bacterium behavior have shown that a simple con- 
nection between the stochastic algorithm of motion and the average chemotactic response 
is far from obvious [9 12] - In particular, it appeared that favorable chemotactic drift could 
not be reconciled with favorable accumulation at long times, and chemotaxis was viewed as 
resulting from a compromise between the two The optimal nature of this compromise 



In various approximations, while the 



in bacterial chemotaxis was examined in Ref. 
negative part of the response kernel was key to favorable accumulation in the steady state, it 
suppressed the drift velocity. Conversely, the positive part of the response kernel enhanced 
the drift velocity but reduced the magnitude of the chemotactic response in the steady state. 

Here, we carry out a detailed study of the chemotactic behavior of a single bacterium in 
one dimension. We find that, for an 'adaptive' response kernel (i.e., when the positive and 
negative parts of the response kernel have equal weight such that the total area under the 
curve vanishes), there is no incompatibility between a strong steady-state chemotaxis and 
a large drift velocity. A strong steady-state chemotaxis occurs when the positive peak of 
the response kernel occurs at a time much smaller than tr and the negative peak at a time 
much larger than tr, in line with experimental observation. Moreover, we obtain that the 
drift velocity is also large in this case. For a general 'non- adaptive' response kernel (i.e., 
when the area under the response kernel curve is non- vanishing) , however, we find that a 
large drift velocity indeed opposes chemotaxis. Our calculations show that, in this case, a 
position-dependent diffusivity is responsible for chemotactic accumulation. 

In order to explain our numerical results, we propose a simple coarse-grained model which 
describes the bacterium as a biased random walker with a drift velocity and diffusivity, both 
of which are, in general, position-dependent. This simple model yields good agreement with 



results of detailed simulations. We emphasize 



grained descriptions of E. coli chemotaxis 14 



rat our model is distinct from existing coarse- 



-Il7|. In these, coarse- graining was performed 



over left- and right-moving bacteria separately, after which the two resulting coarse-grained 
quantities were then added to obtain an equation for the total coarse-grained density. We 
point out why such approaches can fail and discuss the differences between earlier models 
and the present coarse-grained model. 



Models 



Following earlier studies of chemotaxis 10|, llSj, we model the navigational behavior of a 



bacterium by a stochastic law of motion with Poissonian run durations. A switch from run 
to tumble occurs during the small time interval between t and t + dt with a probability 

dt 



{1-T[c}}. 



(1) 
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Here, tr ~ Is and J 7 [c] is a functional of the chemical concentration, c(f), experienced by 
the bacterium at times t' < t. In shallow nutrient gradients, the functional can be written 

as 

7 [c] = j dt'R(t - t')c{t') (2) 

J — oo 

The response kernel, R(t), encodes the action of the biochemical machinery that processes 
input signals from the environment. Measurements of the change in the rotational bias of 
a flagellar motor in wild -typ e bacteria, in response to instantaneous chemoattractant pulses 



were reported in Refs. 



18 



19l | ; experiments were carried out with a tethering assay. The 



response kernel obtained from these measurements has a bimodal shape, with a positive 
peak around t ~ 0.5s and a negative peak around t ~ 1.5s (see Fig. [1]). The negative lobe 
is shallower than the positive one and extends up to t ~ 4s, beyond which it vanishes. The 
total area under the response curve is close to zero. As in other studies of E. coli chemotaxis, 
we take this response kernel to describe the modulation of run duration of swimming bacteria 
j912|. Recent experiments suggest that tumble durations are not modulated by the chemical 



environment and that as long as tumbles last long enough to allow for the reorientation of 
the cell, bacteria can perform chemotaxis successfully 20l. |21|. 

The model defined by Eqs. [T]and|2]is linear. Early experiments pointed to a non-linear, in 



effect a threshold-linear, behavior of a bacterium in response to chemotactic inputs 18l. Il9|. 
In these studies, a bacterium modulated its motion in response to a positive chemoattractant 
gradient, but not to a negative one. In the language of present model, such a threshold-linear 
response entails replacing the functional defined in Eq. [2] by zero whenever the integral is 
negative. More recent experiments suggest a different picture, in which a non-linear response 
is expected only for a strong input signal whereas the response to weak chemoattractant 
gradient is well described by a linear relation [22]. Here, we present an analysis of the linear 
model. For the sake of completeness, in Supporting Information, we present a discussion 
of models which include tumble modulations and a non-linear response kernel. Although 
recent experiments have ruled out the existence of both these effects in E.coli chemotaxis, 
in general such effects can be relevant to other systems with similar forms of the response 
function. 

The shape of the response function hints to a simple mechanism for the bacterium to 
reach regions with high nutrient concentration. The bilobe kernel measures a temporal 
gradient of the nutrient concentration. According to Eq. (TJ if the gradient is positive, runs 
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are extended; if it is negative, runs are unmodulated. However, recent literature 0, Q, 12] 
has pointed out that the connection between this simple picture and a detailed quantitative 
analysis is tenuous. For example, de Gennes used Eqs. H]to calculate the chemotactic drift 
velocity of bacteria [9]. He found that a singular kernel, R(t) = a8(t — A), where S is a Dirac 
function and a a positive constant, lead to a mean velocity in the direction of increasing 
nutrient concentration even when bacteria are memoryless (A = 0). Moreover, any addition 
of a negative contribution to the response kernel, as seen in experiments (see Fig. [1]), lowered 
the drift velocity. Other studies considered the steady-state density profile of bacteria in 
a container with closed walls, both in an approximation in which correlations between run 



durations and probability density were ignored 12| and in a n ap proximation in which the 
memory of the bacterium was reset at run-to-tumble switches Both these studies found 
that, in the steady state, a negative contribution to the response function was mandatory 
for bacteria to accumulate in regions of high nutrient concentration. These results seem to 
imply that the joint requirement of favorable transient drift and steady-state accumulation is 
problematic. The paradox was further complicated by the observation [h| that the steady- 
state single-bacterium probability density was sensitive to the precise shape of the kernel: 
when the negative part of the kernel was located far beyond tr it had little influence on the 



steady-state distribution 12|. In fact, for kernels similar to the experimental one, model 
bacteria accumulated in regions with low nutrient concentration in the steady state flo| . 



Results 



Simulations and analytical treatment of chemotactic bacterial accumulation 

In order to resolve these paradoxes and to better understand the mechanism that leads to 
favorable accumulation of bacteria, we perform careful numerical studies of bacterial motion 
in one dimension. In conformity with experimental observations Q] , we do not make 
any assumption of memory reset at run-to-tumble switches. 

We model a bacterium as a one-dimensional non-Markovian random walker. The walker 
can move either to the left or to the right with a fixed speed, v, or it can tumble at a given 
position before initiating a new run. In the main paper, we present results only for the case 
of instantaneous tumbling with tt = 0, while results for non- vanishing tt are discussed in 
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Supporting Information. There, we verify that for an adaptive response kernel tt does not 
have any effect on the steady-state density profile. For a non-adaptive response kernel, the 
correction in the steady-state slope due to finite tt is small and proportional to r^/r^. 

The run durations are Poissonian and the tumble probability is given by Eq. [TJ The 
probability to change the run direction after a tumble is assumed to have a fixed value, 
q, which we treat as a parameter. The specific choice of the value of q does not affect 
our broad conclusions. We find that, as long as q 7^ 0, only certain detailed quantitative 
aspects of our numerical results depend on q. (See Supporting Information for details on this 
point.) We assume that bacteria are in a box of size L with reflecting walls and that they 
do not interact among each other. We focus on the steady-state behavior of a p opulation. 



Reflecting boundary conditions are a simplification of the actual behavior 23|, |24j ; as long as 



the total 'probability current' (see discussion below) in the steady state vanishes, our results 
remain valid even if the walls are not reflecting. 

As a way to probe chemotactic accumulation, we consider a linear concentration profile 
of nutrient: c(x) = ex. We work in a weak gradient limit, i.e., the value of ac is chosen 
to be sufficiently small to allow for a linear response. Throughout, we use c = 1/L in our 
numerics. From the linearity of the problem, results for a different attractant gradient, k/L, 
can be obtained from our results through a scaling factor k. In the linear reigme, we obtain a 
spatially linear steady-state distribution of individual bacterium positions, or, equivalently, 
a linear density profile of a bacterial population. Its slope, which we denote by /3, is a 
measure of the strength of chemotaxis. A large slope indicates strong bacterial preference 
for regions with higher nutrient concentration. Conversely, a vanishing slope implies that 
bacteria are insensitive to the gradient of nutrient concentration and are equally likely to be 
anywhere along the line. We would like to understand the way in which the slope /3 depends 
on the different time scales present in the system. 



Results with non-adaptive response kernels 

One particular advantage of a linear model is that a general problem can be solved by su- 
perposing the solutions of simpler problems — namely, with delta-function response kernels — 
with suitably chosen coefficients. Thus, solving the problem with a singular response kernel 
amounts to a full solution and we focus here on this case. 



8 



In our simulations, we start from an arbitrary bacterium position inside a box of size L. 
Each time step has a duration dt, during which a running bacterium moves over a distance 
vdt. This distance corresponds to one lattice spacing in our model, in which a lattice is 
introduced because time is discretized. Throughout the numerics, we use dt = 0.01s and 
v = lO/zm/s, which means that the lattice spacing in our simulations is 0.1/zm. Results for 
different values of v can be obtained by rescaling the lattice spacing accordingly. At the 
end of each time step, we compute the functional defined in Eq. [2j for a singular response 
kernel, R(t) = aS(t — A), this takes the form ac [x(t — A)], where c [x(t — A)] is the nutrient 
concentration experienced by the bacterium at time t — A. At the end of each time step 
the bacterium either tumbles, with a probability (1 — ac [x(t — A)]) dt/r^, or continues to 
move in the same direction. At every tumble, the velocity of the bacterium is reversed with 
a probability q. 

The system reaches a steady state over a time scale which is of order L 2 /D, where the 
diffusivity is given by D = v 2 tr. We verify numerically that after this time the bacterial 
density profile inside the box does not change further and assumes a time-independent 
linear form. We focus on the slope, 0, of this profile. For an experimental realization of 
the steady-state behavior of a single bacterium, we provide here an estimate of the time 
scales and length scales involved. Since the long-time behavior of the system is diffusive 
(see the discussion of the coarse-grained model below), the relaxation time is L 2 /D. Our 
results on the steady-state distribution of bacteria hold, realistically, if this relaxation time 
does not exceed the typical division time of an E. coli bacterium, which is of the order of 
30 minutes. Substituting experimental values for the parameters, we find the description 
should be valid for system sizes L < 400/im. In our simulations, we use a somewhat larger 
system (L = lOOOyum) so as to have cleaner results with negligible effects of the reflecting 
walls at the two boundaries. (Numerics data show that the width of the boundary layer is 
about ~ 80/im.) 

According to our numerical simulations, for a < 0, /3 increases with A and displays a 
plateau for A ^> tr (Fig. [2). Simulations probing various values of tr also confirmed that 
(3 = F(A/tr), i.e., that the slope is a scaling function of A/tr. Clearly, for positive a the 



sign of f3 is simply reversed, which corresponds to an unfavorable chemotaxis [12l. Il5|. 

For small A, one can write down an approximate master equation for left-mover and right- 
mover densities and use it to show that the slope increases linearly with A (see Supporting 
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FIG. 2: The slope /3 (scaled by a factor of 10 8 ) as a function of A/77?, for the choice of response 
kernel R(t) = aS(t — A). Note that for A ^> tr the slope saturates to a non- vanishing value. The 
symbols +, *, and □ correspond to slopes measured in simulations with tr = 1.5, 1.0, and 0.5 
seconds, respectively. The black solid circles are derived from our coarse-grained formulation (Eq. 
E]). Here q = 0.5, a = -0.02, L = 1000/im, c = 0.001/im" 1 , v = W^m/s. 

Information for details). It is surprising, however, that the slope appears to saturate to 
a non-vanishing value for A ^> tr. Indeed one would expect that, if the response kernel 
relies on a time much earlier than t — tr, a large enough number of tumbles occur between 
this past time and the present time so as to eliminate any correlation between the nutrient 
concentration in the past and the present direction of motion. If this argument holds, one 
would expect that the slope /3 vanish for A ^> tr. Below, we return to this argument and 
explain why it is misleading. 

Results with adaptive response kernels 




For wild-type bacteria, the total area under the response kernel vanishes (Fig. [T]). As a 
result, their behavior is adaptive: chemotaxis is insensitive to the overall level of nutrient, 
but sensitive to spatial variations 18|, [l9| . In this section, before examining the case of a 
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bilobe response kernel similar to the experimental one, we consider a toy model defined by 
the difference of two singular forms: R(t) = aS(t — Ai) — aS(t — A 2 ), with a > 0. Because 
our problem is linear, the steady-state slope of bacterial density, (3, can be calculated from 
a simple linear superposition, as: 



Since the function F(-) is monotonic, the absolute value of /3 increases with the difference 
of Ai and A 2 . Strong chemotaxis occurs when Aj = and A 2 3> tr. 

We now turn to the experimental case of a bilobe response kernel. It is not computa- 
tionally feasible to work with the complete set of experimental data [19|, so we have used a 
discrete subset (Fig. [1]) which we represent as a series of delta-functions. Given this approx- 
imate response kernel, we investigate the behavior of the slope as a function of tr. Based 
on our results for the case of two delta functions, we expect that chemotaxis be weak if tr is 
either much smaller than the delay of the positive peak in the response kernel or much larger 
than the delay of the negative peak. We expect optimum chemotaxis for a value of tr that 
falls in between the two delays. We verify this prediction in Fig. [3] (in the linear model). We 
note that the maximum slope occurs for a value of tr close to the experimentally recorded 
value of about Is. 



In order to gain insight into our numerical results, we developed a simple coarse-grained 
model of chemotaxis. For the sake of simplicity, we first present the model for a non-adaptive, 
singular response kernel, R(t) = aS(t — A), and, subsequently, we generalize the model to 
adaptive response kernels by making use of linear superposition. 

The memory trace embodied by the response kernel induces temporal correlations in 
the trajectory of the bacterium. However, if we consider the coarse-grained motion of the 
bacterium over a spatial scale that exceeds the typical run stretch and a temporal scale 
that exceeds the typical run duration, then we can assume that it behaves as a Markovian 
random walker with drift velocity V and diffusivity D. Since the steady-state probability 




(3) 



Coarse-grained description of chemotaxis as diffusion with drift 
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FIG. 3: The slope /3 (scaled by a factor of 10 8 ) as a function of tr for the experimental response 
kernel shown in Fig. [TJ Open squares: numerical results from simulations. Solid circles: prediction 
of the coarse-grained model. Here, q = 0.4, L = lOOO^m, c = O.OOl^m" 1 , v = 10/um/s. 

distribution, P(x) = P(A, tr, x), is flat for a = 0, for small a we can write 

P = P + aV(A, tr, x) + o(a 2 ), (4) 
D = D + aV(A, tr, x) + o(a 2 ), (5) 
V = aV{A,T R ,x) + o(a 2 ). (6) 

Here, Pq = 1/L and D = v 2 tr. Since we are neglecting all higher order corrections in a, 
our analysis is valid only when a is sufficiently small. In particular, even when A 3> tr, we 
assume that the inequality A/tr 1/at is still satisfied. The chemotactic drift velocity, V, 
vanishes if a = 0; it is defined as the mean displacement per unit time of a bacterium starting 
a new run at a given location. Clearly, even in the steady state when the current J, defined 
through dfP = —d x J, vanishes, V may be non- vanishing (see Eq. [S] below). In general, 
the non-Markovian dynamics make V dependent on the initial conditions. However, in the 
steady state this dependence is lost and V can be calculated, for example, by performing 
a weighted average over the probability of histories of a bacterium. This is the quantity 
that is of interest to us. An earlier calculation by de Gennes showed that, if the memory 
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preceding the last tumble is ignored, then for a linear profile of nutrient concentration 
the drift velocity is independent of position and takes the form V = acv 2 tr exp(— A/tr) 
While the calculation applies strictly in a regime with A < tj (because of memory 
erasure), in fact its result captures the behavior well over a wide range of parameters (see 
Fig. Hj). To measure V in our simulations, we compute the average displacement of the 
bacterium between two successive tumbles in the steady state, and we extract therefrom the 
drift velocity. (For details of the derivation, see Supporting Information.) We find that V 
is negative for a < and that its magnitude falls off with increasing values of A (Fig. H]). 
We also verify that V indeed does not show any spatial dependence (data shown in Fig. 
S3 of Supporting Information). We recall that, in our numerical analysis, we have used a 
small value of a; this results in a low value of V. We show below that for an experimentally 
measured bilobe response kernel, obtained by superposition of singular response kernels, the 
magnitude of V becomes larger and comparable with experimental values. 
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FIG. 4: The chemotactic drift velocity, V, as a function of A, for the response kernel R(t) = 

n 

a5(t — A). Solid circles: numerical results. Line: approximate analytical results from [9]. tr = Is 
and other numerical parameters as in Fig. [2] 

To obtain the diffusivity, D, we first calculate the effective mean free path in the coarse- 
grained model. The tumbling frequency of a bacterium is (1 — acx(t — A))/tr and depends 



T 
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on the details of its past trajectory. In the coarse-grained model, we replace the quantity 
acx(t — A) by an average ac(x(t — A)) over all the trajectories within the spatial resolution 
of the coarse-graining. Equivalently, in a population of non-interacting bacteria, the average 
is taken over all the bacteria contained inside a blob, and, hence, {x(t — A)) denotes the 
position of the center of mass of the blob at a time t — A in the past. As mentioned 
above, the drift velocity is proportional to a, so that ac(x(t - A)) = acx(t) + 0(a 2 ). The 
average tumbling frequency then becomes (1 — acx))/r R and, consequently, the mean free 
path becomes t mfp = r R /(l — acx) ~ t r (1 + acx). As a result, the diffusivity is expressed 
as D = v 2 tmfp — v 2 tr{1 + acx). We checked this form against our numerical results (Fig. 
E]). 

Having evaluated the drift velocity, V, and the diffusivity, D, we now proceed to write 



down the continuity equation (for a more rigorous but less intuitive approach, see 
For a biased random walker on a lattice, with position-dependent hopping rates d + (x) and 
d~(x) towards the right and the left, respectively, one has V = a (d + (x) — d~(x)) and 
D = a 2 (d + (x) + d~(x)) /2, where a is the lattice constant. In the continuum limit, the 
temporal evolution of the probability density is given by a probability current, as 



where the current takes the form 



d t P = -d x J, (7) 



J = VP-d x (DP) . (8) 



For reflecting boundary condition, J = in the steady state. This constraint yields a 
steady-state slope 



= ad x P = a^(V- d x V) = -r-f- - — (9) 



:— V - d x V) = — - 

V Lv z t r L 

for small a. We use our measured values for V and D (Figs. H]and[5]), and compute the 
slope using Eq. [9j (For details of the measurement of V, see Supporting Information.) We 
compare our analytical and numerical results in Fig. [21 which exhibits close agreement. 

According to Eq. [HJ steady-state chemotaxis results from a competition between drift 
motion and diffusion. For a < 0, the drift motion is directed toward regions with a lower 
nutrient concentration and hence opposes chemotaxis. Diffusion is spatially dependent and 
becomes small for large nutrient concentrations (again for a < 0), thus increasing the ef- 
fective residence time of the bacteria in favorable regions. For large values of A, the drift 



14 



E 

Q 

>* 

> 
en 



99.2 
99 
98.8 
98.6 
98.4 
98.2 
98 
97.8 




200 300 400 500 600 700 800 

position x (urn) 



FIG. 5: The diffusivity, D(x), as a function of position, x, for the response kernel R(t) = a5(t — A) 
with A = Is (+) and 2s (x). Instead of plotting D{x) for the entire range of x, we leave out 
boundary regions to avoid the effect of the reflecting walls. (From the numerics, the width of 
the boundary layer is ~ 80/um.) D(x) falls off linearly with x and is independent of A. Data 
fitting yields D(x) = 99.5 — 0.00197x and the coarse-grained model predicts D(x) = v 2 tr(1 + acx). 
For the chosen set of parameters, v 2 tr = 100/xto 2 /s and the v 2 TRac = —0.002. The discrepancy 
between the numerical and the predicted slopes is due to higher-order corrections in a, while 
discretization of space in simulations causes the slight mismatch in the constant term, tr = Is 
and other numerical parameters are as in Fig. [21 

velocity vanishes and one has a strong chemotaxis as A increases (Fig. [2]). Finally, for 
A = 0, the calculation by de Gennes yields V = acv 2 T^ which exactly cancels the spatial 
gradient of D (to linear order in a), and there is no accumulation [9], 

These conclusions are easily generalized to adaptive response functions. For R(t) = 
a5(t — A x ) — a5(t — A 2 ), within the linear response regime, the effective drift velocity 
and diffusivity can be constructed by simple linear superposition: The drift velocity reads 
V = ctV^Ax) — aV(A 2 ). Interestingly, the spatial dependence of D cancels out and D = 
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Dq = v 2 r R . The resulting slope then depends on the drift only and is calculated as 



In this case, the coarse-grained model is a simple biased random walker with constant 
diffusivity. For Ai < A 2 and a > 0, the net velocity, proportional to a(V(Ai) — V(A 2 )), 
is positive and gives rise to a favorable chemotactic response, according to which bacteria 
accumulate in regions with high food concentration. Moreover, the slope increases as the 
separation between Ax and A 2 grows. We emphasize that there is no incompatibility between 
strong steady-state chemotaxis and large drift velocity. In fact, in the case of an adaptive 
response function, strong chemotaxis occurs only when the drift velocity is large. 

For a bilobe response kernel, approximated by a superposition of many delta functions 
(Fig. CP, the slope, (3, can be calculated similarly and in Fig. [3] we compare our calculation 
to the simulation results. We find close agreement in the case of a linear model with a 
bilobe response kernel and, in fact, also in the case of a non-linear model (see Supporting 
Information). 

The experimental bilobe response kernel R(t) is a smooth function, rather than a finite 
sum of singular kernels over a set of discrete A values (as in Fig. [1]). Formally, we integrate 
singular kernels over a continuous range of A to obtain a smooth response kernel. If we 
then integrate the expression for the drift velocity obtained by de Gennes, according to this 
procedure, we find an overall drift velocity V ~ 0.3 fim/s, for the concentration gradient 
considered (Vc = O.OOl/im -1 ). By scaling up the concentration gradient by a factor of k, 
the value of V can also be scaled up by k and can easily account for the experimentally 
measured velocity range. 



We carried out a detailed analysis of steady-state bacterial chemotaxis in one dimension. 
The chemotactic performance in the case of a linear concentration profile of the chemoat- 
tractant, c(x) = cx, was measured as the slope of the bacterium probability density profile 
in the steady state. For a singular impulse response kernel, R(t) = a5(t — A), the slope was 
a scaling function of A/tr, which vanished at the origin, increased monotonically, and sat- 
urated at large argument. To understand these results we proposed a simple coarse-grained 




(10) 
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model in which bacterial motion was described as a biased random walk with drift velocity, 
V, and diffusivity, D. We found that for small enough values of a, D was independent of A 
and varied linearly with nutrient concentration. By contrast, V was spatially uniform and 
its value decreased monotonically with A and vanished for A ^> tr. We presented a simple 
formula for the steady-state slope in terms of V and D. The prediction of our coarse-grained 
model agreed closely with our numerical results. Our description is valid when a is small 
enough, and all our results are derived to linear order in a. We assume A/tr <C I /a is 
always satisfied. 

Our results for an impulse response kernel can be easily generalized to the case of response 
kernels with arbitrary shapes in the linear model. For an adaptive response kernel, the spatial 
dependence of the diffusivity, D, cancels out but a positive drift velocity, V, ensures bacterial 
accumulation in regions with high nutrient concentration, in the steady state. In this case, 
the slope is directly proportional to the drift velocity. As the delay between the positive and 
negative peaks of the response kernel grows, the velocity increases, with consequent stronger 
chemotaxis. 

Earlier studies of chemotaxis 
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171 ] put forth a coarse-grained model different from ours. 
In the model first proposed by Schnitzer for a single chemotactic bacterium jla |. he argued 
that, in order to obtain favorable bacterial accumulation, tumbling rate and ballistic speed 
of a bacterium must both depend on the direction of its motion. In his case, the continuity 
equation reads 



d t P = d x 



IlVr - IrVl p 2 Vr + Vl d ( VRVL P 



(11) 



1L+1R 1R + 1L \Vr + V L 

where vl(r) is the ballistic speed and 'Jl(r) is the tumbling frequency of a bacterium moving 
toward the left (right). For E. coli, as discussed above, Vl = vr = v, a constant independent 
of the location. In that case, Eq. [TT] predicts that in order to have a chemotactic response 
in the steady state, one must have a non-vanishing drift velocity, i.e., (■Jl v r — IrVl) /(jl + 
7k) 7^ 0. This contradicts our findings for non- adaptive response kernels, according to 
which a drift velocity only hinders the chemotactic response. The spatial variation of the 
diffusivity, instead, causes the chemotactic accumulation. This is not captured by Eq. [UJ 
In the case of adaptive response kernels, the diffusivity becomes uniform while the drift 
velocity is positive, favoring chemotaxis. Comparing the expression of the flux, J, obtained 
from Eqs. [7] and [8] with that from Eq. [Til an d matching the respective coefficients of P 
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and d x P, we find D = 2v R v L /('j R + j L ) and V = (j L v R - j R v L )/(j L + j R ). As we argued 
above in discussing the coarse-grained model for adaptive response kernels, both D and V 
are spatially independent. This puts strict restrictions on the spatial dependence of vl(r) 
and 7l(_r). For example, as in E. coli chemotaxis vl = v R = v, our coarse-grained description 
is recovered only if jl and j R are also independent of x. 

We comment on a possible origin of the discrepancy between our work and earlier treat- 



ments. In Ref. 15], a continuity equation was derived for the coarse-grained probability 
density of a bacterium, starting from a pair of approximate master equations for the prob- 
ability density of a right-mover and a left-mover, respectively. As the original process is 
non-Markovian, one can expect a master equation approach to be valid only at scales that 
exceed the scale over which spatiotemporal correlations in the behavior of the bacterium 
are significant. In particular, a biased diffusion model can be viewed as legitimate only if 
the (coarse-grained) temporal resolution allows for multiple runs and tumbles. If so, at the 
resolution of the coarse-grained model, left- and right-movers become entangled, and it is 
not possible to perform a coarse-graining procedure on the two species separately. Thus one 
cannot define probability densities for a left- and a right-mover that evolves in a Markovian 
fashion. In our case, left- and right-movers are coarse-grained simultaneously, and the total 
probability density is Markovian. Thus, our diffusion model differs from that of Ref. [ijj] 
because it results from a different coarse-graining procedure. The model proposed in Ref. 



3 



has been used extensively to investigate collective behaviors of E. coli bacteria such as 



pattern formation Il71 |. It would be worth asking whether the new coarse-grained 



description can shed new light on bacterial collective behavior. 
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Supporting Information 

Small- argument behavior of F (A/tr) 

Here, we argue that for small A/tr the function F(A/tr) is linear. First, note that for 
an impulse response kernel, R(t) = a5(t — A) the tumbling probability during the interval 
[t, t + dt] is [1 — acxit — A)] cH/tr. For small A, we can assume that no tumbling occurs 
during the interval [t — A, t]. Then the effective tumbling rates become [1 — ac(x — vA)] /tr 
for right-movers and [1 — ac(x + v A)] /tr for left-movers. Based on this observation, we can 
write a pair of master equations that govern the densities of left-movers and right-movers, 
and we can solve them in the steady state. We obtain the slope, (3, expressed as 

*=-^. (12) 
Thus, in the limit A <C tr the scaling function F(A/tr) depends linearly on its argument. 



Chemotactic drift velocity as a function of position 

Here, we show a sample of our numerical results on the drift velocity V as a function 
of position. We measure the average displacement of a bacterium in the steady state (see 
the corresponding discussion in the paper), and we compute V therefrom. Specifically, we 
perform the following procedure. We denote by m^x) and m^(x) the total number of 
leftward and rightward runs, respectively, that initiate at the position x, at a time t, for a 
population of non-interacting bacteria. These quantities are well-defined in the context of 
our numerics because space is discretized. Furthermore, let Sl(x) and Sr(x) be the total 
leftward and rightward displacement, respectively, of the bacteria that undergo these runs. 
The average displacement per run is then given by [Sr(x) — SL(x)]/[mR(x) + m^x)]. Of 
course, the choice of the time scale as the duration of a run is arbitrary and other choices are 
equally valid. To leading order, this average displacement is linear in a. The drift velocity 
V, to order a, is then obtained by dividing this average displacement by tr. (Note that any 
0(a) correction to tr leads to 0(a 2 ) correction to V, which we ignore here.) We find that 
up to the noise present in our numerical measurement, V does not display any dependence 
on position. An example is illustrated in Fig. [61 
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FIG. 6: Drift velocity, V, as a function of position, x, for the case of a singular response kernel 
R{t) = aS(t — A) (red line). Instead of showing the plot for the entire range of x, we leave out 
boundary regions to discard the effect of the reflecting walls. Our numerics show that the width of 
the boundary layer is ~ 80/im. Here, L = 1000/im, q = 1.0, tr = Is, A = Is, tt = 0, v = Wfim/s, 
a = —0.1, c = O.OOlfim^ 1 . Based on the data shown, the drift velocity is 0.001 ±0.0001/im/s (blue 
line). 

Dependence of V on the turning probability q 

In our model, q denotes the turning probability, i.e., the probability that the run direction 
inverses after a tumble. Our numerical explorations indicate that changing the value of q 
does not affect the qualitative behavior of the system. However, the numerical value of the 
drift velocity, and the value of /3 depend on q. In Fig. [TJ we exhibit the variation of V as a 
function of q, in the case of an adaptive response function. 

Results in the non-linear model 

Some earlier experiments indicate that bacteria modulate their run durations in response 
to a positive concentration gradient, but not to a negative one. In order to incorporate this 
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FIG. 7: The drift velocity, V , as a function of the turning probability, q, for the case of a singular 
adaptive response kernel R(t) = a5(t — Ai) — ad(t — A2). Here, L = WOOfim, tr = Is, Ai = 0.5s, 
A2 = 1.5s, tt = 0, v = Wfim/s, a = 0.1, c = O.OOl/im -1 . 

feature in our model, we have to go beyond the linear response regime. In the non-linear 
model, whenever the linear functional (Eq. 2) becomes negative, it is replaced by 0. This 
is the only difference with the linear model. Thus, for a purely positive response kernel the 
non-linear model behaves identically to the linear model, while for a purely negative response 
kernel the non-linear model displays no chemotaxis whatsoever. Hereafter, we examine only 
adaptive response kernels with balanced positive and negative contributions. 

We first consider the idealized response kernel made of the superposition of positive and 
negative delta functions, R(t) = aS(t — Ai) — a8(t — A2). Simulations show that many 
qualitative features of the linear model still hold in the non-linear model. The scaling form 
valid in the linear case breaks down in the non-linear case, but the slope /3 increases with the 
separation between Ai and A2 and ultimately saturates to a non-vanishing value. Figures 
|8] and M display results of simulations. As expected, strong chemotaxis occurs when Ai = 
and A 2 is substantially larger than tr. We have also verified that tumbling does not have 
much of an effect on the slope. 

In the experimental case of a bilobe response kernel (Fig. 1 in the paper), we find that 



_j 1 1 1 1 1 1 i_ 



23 



2.6 | , , , , , , , 1 

2.4 - • - 

2.2 - 

2 - 

I 1.8 - 

CO. 

1.6 - 
1.4 - 
1.2 - 

1 ' ' ' ' ' ' ' ' 

0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 

A 2 (S) 

FIG. 8: The slope, /3 (scaled by a factor of 10 8 ), as a function of A2, for Ai = 0.5s, in a non-linear 
model with balanced response kernel, R(t) = a5{t — Ai) — aS(t — A2). As in the linear model, 
/3 increases with the difference of Ai and A2. Here, tr = 1 sec, q = 0.4, a = 0.1, L = 1000/im, 
c = 0.001/xm -1 , v = Wfim/s. 

strong chemotaxis occurs when tr lies between the positive and the negative peaks of the 
response kernel, as found in the linear case. For smaller or larger values of tr, chemotaxis 
becomes weak. Figures [TD] and [TT] show our numerical results for q = 0.4 and q = 0.5, 
respectively. We note that, in both plots the value of tr for which the slope is maximum 
falls close to the experimental value of about Is. However, the exact position of the maximum 
depends on q. For q = 0.4 (Fig JTOj) maximum occurs at tr — 0.8s, while for q = 0.5 (Fig. 
[TTT) the maximum occurs at tr c± Is. 

Chemotaxis with non-vanishing tumbling durations 

During a tumbling event the bacterium rotates about itself in a random fashion without 
any significant displacement. In a homogeneous nutrient concentration the average tumbling 
duration is 0.1s, which is much smaller than the average run duration of Is. In the steady 
state one therefore expects that the bacterium spends only a fraction t^/tr ~ 0.1 of the time 
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FIG. 9: The slope, /3 (scaled by a factor of 10 s ), as a function of Ai, for A2 = 1.5s, in a non-linear 
model with balanced response kernel, R(t) = ad(t — Ai) — a5(t — A2). Numerical parameters are 
as in Fig. [HJ 

in the tumbling state. For this reason, studies of chemotaxis often assume instantaneous 
tumbling. 

It was shown recently that the existence of non-vanishing tumbling duration can yield 
interesting results: even with a punctual response kernel, R(t) = aS(t — A) with A = 0, i.e. 
a memoryless bacterium, one can observe a chemotactic response (Kafri et ai., 2008). Here, 
we provide a simplified derivation of the steady-state density profile in this Markovian limit. 
The result will prove useful also for the analysis of the non-Markovian case with A ^ 0. 

Let L(x,t) and R(x,t) be the density of left-movers and right-movers, respectively, at 
location x and time t. We denote by T R (x,t) and T L (x,t) the densities of tumblers that 
were moving to the right and left, respectively, before tumbling. For A = 0, the time 
evolution of these quantities can be described by master equations. In the case in which 
tumble durations are not modulated and tumble-to-run switches always occur at a fixed rate 
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FIG. 10: The slope, /3 (scaled by a factor of 10 s ), as a function of tr, in a non-linear model with 
the experimental bilobe response kernel of Fig. 1 in the paper. Open squares: numerical results 
from simulations. Solid circles: prediction of the coarse-grained model. Numerical parameters as 
in Fig. 3 in the paper. The inset shows the steady-state density profiles of the bacterial population 
for tr = 0.3,0.8s (red and blue curve), respectively. 

1/tt, the master equations read 

d ^A= - vdm R{x,t)+T*(x,t)^ 

+T L (x,t)±-R(x,ty-^, (13) 

dJ ^= vd x L(x,t)+T^x,t)^ 

+T R (x,t)±-L(x,t) 1 -^, (14) 

R(x,t)l=*2-T*(x,t)± (15) 



dT R (x,t) 
dt 

dT L (x,t) 



L(x,ty-^-T L (x,t)±. (16) 



{ft T Q T~T 

We consider perfectly reflecting boundary conditions at x = and x = L. This implies that, 
in the steady state, we must have R(x) = L(x). The steady-state (total) density at location 
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FIG. 11: The slope, /3 (scaled by a factor of 10 8 ), as a function of tr, in a non-linear model with the 
experimental bilobe response kernel of Fig. 1 in the paper. Here, q = 0.5 and the other numerical 
parameters are as in Fig. 3 in the paper. Comparison with Fig. [10] shows that the position of 
the maximum depends on the value of q. Open squares: numerical results from simulations. Solid 
circles: prediction of the coarse-grained model. 



N(x)= R{x)+L{x)+T R {x) + T L \ 



x 



x then becomes 



= i^k^^ + ^i 1 -^) ■ ( 17 ) 

Therefore, the slope of the steady-state density profile is given by 

(3 = G (—} = -2ac\ T -^—± ( (18) 

VW Lr 2 + ^(2-a) 1 } 

In Fig. [12] we compare this result with the slope measured in simulations. This plot 
demonstrates that, even in the absence of any memory or modulation of tumbling durations, 
it is possible to obtain chemotaxis in the steady state. 

This result for the slope is slightly different for the case in which tumbling durations are 
modulated: then, tt in the above master equations is replaced by t t /{1 + acx). Solving for 
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FIG. 12: The steady-state slope (3 as a function o/tt/to for A = and unmodulated tumbling. The 
solid line corresponds to the exact result from Eq. Ilffl Here we have used L = lOOO^m, a = —0.1, 
q = 0.4, c = 0.001/um -1 , v = 10 [im/ sec. 

the steady state, we find a total density 



N(x) = R(x) + L{x) + T R {x) + T L \ 

i l—acx 



X 



(19) 



For acx 1, this is approximated by a linear form and the slope becomes 

t t \ Icxctt 1 



= G[!Z)=- ^,* „ v (20) 



TO 



We compare this analytical result with simulations in Fig. [131 which shows a systematic de- 
viation for large argument. We have verified that this mismatch originates in the linearizing 
approximation step from Eq. [19] to Eq. [201 

For A ^ and when both runs and tumbles are modulated, we measure the density 
profile in numerical simulations. For R(t) = a5(t — A), our numerics indicate that the 
steady-state slope, (3, is a sum of the Markovian component, G, defined in Eq. [201 an d a 
non-Markovian component, F, which depends on A/t# but is independent of tt'- 

P = f(^-) + g( T A. (21) 



trJ \ t r 
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FIG. 13: The slope, f3, as a function of tt/tr for A = and modulated tumbling. Simulation 
parameters are L = 1000/im, a = —0.1, q = 0.4, c = 0.001/im -1 , v = 10//ra/s. 

In Fig. [TU we exhibit this scaling form as a function of t^/tr for a fixed, non- vanishing 
value of A/tr. Our results suggests that (3 is made up of two contributions: one from the 
modulating runs, encoded in F, and one from non- instantaneous tumbles, encoded in G. 
The latter contribution is independent of A. 

Finally, we can infer more general results from the simple form in Eq. |2TJ In particular, 
for an adaptive response function in the linear model, the positive and negative parts of the 
response function cancel out the effect of non-vanishing tumble durations. In this case, the 
steady-state slope, (3, becomes independent of tt- Interestingly, we find the same applies 
even in the non-linear model (see Fig. H%|) . 
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FIG. 14: The scaling collapse of the slope, f3, as a function of tt / Tj for fixed value o/A/tr. We 
have used A/77? = 0.5 here. The other simulation parameters are as in Fig. \1SX We also plot the 
slope, j3, in the non-linear model with R(t) = a [S(t — Ai) — S(t — A2)], Ai = 0s and A2 = 1.5s 
and tr = Is: f3 does not show any significant dependence on tt- 



